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In this letter, we develop a Newtonian approach to the collapse of galaxy fluctuations of scalar field 
dark matter under initial conditions inferred from simple assumptions. The full relativistic system, 
the so called Einstein-Klein-Gordon, is reduced to the Schrodinger-Newton one in the weak field 
limit. The scaling symmetries of the SN equations are exploited to track the non-linear collapse of 
single scalar matter fluctuations. The results can be applied to both real and complex scalar fields. 
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The most recent cosmological observations pj support 
the Cold Dark Matter (CDM) as the standard model 
of cosmological structure formation. So far, the CDM 
model agrees reasonably well with many observations at 
large and galactic scales |2J, |3|, though its predictions 
at subgalactic scales are still a matter of intense debate 
Pf. However, the nature of CDM remains one of the 
most intriguing problems in modern cosmology. It is 
amazing that CDM can be so predictive and survive the 
confrontation with cosmological data just by assuming 
general features of dark matter such like, for instance, 
that it is made of weakly interacting massive particles 
(WIMPS). On the other hand, with CDM becoming a 
reliable paradigm of structure formation and its nature 
still uncertain, it is possible to test models of dark matter 
as alternatives to the WIMP hypothesis; models in which 
the physical processes can be much more under control. 

One of such alternative models is the scalar field dark 
matter (SFDM), in which dark matter should be com- 
prised of ultra-light scalar field particles [E IE IE IE IE 
Hot ITU . On the one hand, the cosmological evolution 
of SFDM and its linear fluctuations can match those of 
CDM. On the other hand, self-gravitating scalar configu- 
rations can reproduce some general properties of galaxy 
halos as observed today. Going further, it is possible to 
study the non-linear evolution of SFDM fluctuations to 
form gravitationally bound objects, as it was first shown 
in . This is done by numerically evolving the coupled 
Einstein-Klein-Gord on ( EKG) equations for scalar field 
configurations |9l ll2lll3l ll4 | . 

In this letter, our main aim is to study the collapse 
of scalar field dark matter under realistic conditions. In 
this sense, we will try to complete the picture of struc- 
ture formation with scalar field dark matter, from the 
growing of primordial fluctuations up to its non-linear 
collapse. Also, our intention is to draw a guideline as 
simple as possible, for which we will only take into ac- 
count the simplest assumptions; a more detailed picture 
will be presented elsewhere. 

First of all, we consider that SFDM is comprised of 
a minimally coupled real scalar field $ endowed with a 



scalar field potential V(Q), which accounts for the self- 
interactions of the scalar field. Even though the general 
form of the scalar potential can be complicated, we are 
interested only in those potentials that behave as V($>) — 
(ra 2 /2)<I> 2 at late times [a lift LA IS ■ Such a quadratic 
scalar potential is preferred because, in a homogeneous 
and isotropic Universe, the scalar field oscillates coher- 
ently around the minimum of the potential and then the 
scalar energy density p$ ~ m 2 $ 2 scales much the same as 
CDM. Moreover, the scalar field fluctuations in the linear 
regime also grow in the same manner as CDM fluctua- 
tions. Most of the works suggest that a realistic model 
of structure formation should consider a very light mass 
for the boson particles [E IE IE IE IE 03 • We will use here 
the value m = l(T 23 eV = (2.09 yr)" 1 = (0.64 pc)" 1 (we 
are taking units such that c = h = 1 ; also G = m P1 with 
topi is the Planck mass) , which is also a natural scale for 
time and distance within SFDM. 

Due to the similarities of scalar field dark matter with 
the CDM model, it is reasonable to take as a guideline for 
the non- linear regime of structure formation the standard 
spherical collapse model 15]. According to this, a spher- 
ical overdense fluctuation slowly separates out from the 
cosmological expansion, reaches a maximum expansion 
and then collapses under the influence of its own gravity 
only. The important instant here is that of maximum 
expansion, which is called the time of turnaround. 

Before turnaround, we do not expect violent processes 
occurring to the matter fluctuation, and that is why we 
consider the spherical collapse model as a good approx- 
imation for scalar fields. However, we indeed expect 
the gravitational collapse of the matter fluctuation to be 
determined by the intrinsic nature of the scalar fields. 
Hence, the conditions at turnaround will be our initial 
conditions for the gravitational collapse of scalar field 
dark matter. 

According to the spherical collapse model, the non- 
linear density contrast at turnaround is S^ Tn = 
[dp$/p$]( tUTn ^ = 4.55. As we said before, the homoge- 
neous scalar energy density evolves as the standard cold 
dark matter does, and then the local value of the scalar 
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energy density at turnaround is approximately given by 
87rG$ t 2 urn ~ 13.55^0 ,CDM (1 + z turn ) 3 H 2 /m 2 , (1) 



1 is the cur- 
um ~ few is 



where rio,CDM — 0.25, Ho = 70kms _1 Mpc 
rent value of the Hubble parameter, and zt 
the redshift at turnaround. Thus, the local value of the 
scalar field at turnaround is quite small, of the order 
y^<° r c n al) ^5-06><10- 10 . 

With this in mind, we proceed now to describe the 
gravitational colla pse of a fluctuation. As it has been 
shown in [gl [T3 . fl3Lll4|| . the problem is well defined within 
General Relativity, of which we give here a brief descrip- 
tion. Assuming spherical symmetry, the metric is written 
in the form 



ds 2 



a 2 dt 2 + a 2 dr 2 + r 2 (d9 2 + sin 2 9d<p 2 ) , (2) 



where a(t,r),a(t,r) are functions determined by the 
matter distribution through the coupled Einstein-Klein- 
Gordon (EKG) equations, 



G i1u =&-kGT i1v (<$>), D$ = m 2 $. 



(3) 



Here, G^ v is the Einstein tensor corresponding to the 
metric l|2jl. and T^(<&) is the scalar energy- momentum 
tensor. The KG equation arises from the Bianchi identi- 
ties, and □$ = -^d^/^gg^d^]. 

In principle, we should just numerically evolve Eqs. 
(PJ, but as we have to deal with weak gravity we find 
more appropriate to evolve the weak field limit of such 
system of equations, which by the way provides impor- 
tant technical simplifications that let us to have more 
physical insight. In our case, that limit arises when 
a 2 - I, a 2 - 1, V8^G$ < 10" 3 [HQl. 

As it was shown in fllj , such limit is found for complex 
scalar fields through a standard post-Newtonian treat- 
ment, and we now show how a similar procedure can also 
be applied to real scalar fields. We begin by writing such 
real scalar field and the metric coefficients in terms of the 
Newtonian fields ip, U, U2, A, A 2 as 



V8^G$ = e- iT ip(T,x) +C.C, (4) 
a 2 = l + 2U(T,x) + e- 2iT U 2 (T,x) + C.C., (5) 
a 2 = l + 2A(T,x) + e- 2iT A 2 (T,x) + C.C, (6) 

where we have also introduced the dimensionless quan- 
tities t = mi, x = mr. Notice that only U, A are real 
fields. Next, we assume that all the new fields are of order 
C(e 2 ) <C 1, and that their derivatives obey the standard 
post-Newtonian rules d T ~ ed x ~ 0(e 4 ). Therefore, con- 
sidering the leading order terms only, the EKG equations 
now read 



id T ip 

d 2 x (xU) 
d x U 2 



~—d 2 x (x^) + UiP : 
xipip* , 

—x4> 2 ■ 



(7) 

(8) 
(9) 



In addition, A(t,x) = xd x U and A 2 ~ (^(e 4 ), that is, 
the metric coefficient g rr can be taken plainly as time- 
independent. 

That the complete system {7J), |(SJ| and @ is really 
the weak- field manifestation of Eqs. is easily ver- 
ified since the former admits stationary solutions for 
ijj = e~ luJT 4>{x). Under appropriate boundary conditions, 
such solutions are in turn the so-called Newtonian oscilla- 
tons ^(| . Indeed, Eq. 10 only arises for real scalar fields 
and represents the particular oscillatory behavior of the 
metric for oscillatons [Hi. ITiL . 

It should be stressed here that the whole dynamics of 
the EKG system is contained in Eqs. J7J and (JSJ), which 
are the so-called Schrodinger-Newton (SN) equations 0, 
IE 13 ^1 113 j which also stands for the post-Newtonian 
expansion with complex scalar fields 12]. 

Before proceeding further, we need to determine the 
initial scalar field profile. As we said before, prior to 
turnaround, the scalar profile is not determined by its 
own gravity only. In other words, the gravitational po- 
tential U(T tmni x) is not determined by Eq. l|HJl_10]- But, 
-0 should still obey the Schrodinger equation J7J, i.e., the 
initial condition should be in accord with the scalar na- 
ture of the matter fluctuation. 

Thus, the initial condition problem can be solved 
by finding the initial gravitational potential U(0, r) — 
U(tt U rmr). The simplest assumption we can make is to 
approximate the gravitational potential by a spherical 
square well of the form 01 

U(0,r<R ) = U , U(0,r>R ) = 0, (10) 

where Uq = —GMo/Rq is the depth of the gravitational 
well and i?o is taken as the radius containing the total 
initial mass Mq. Next, we look for stationary solutions 
of Eq. Q) in the form -0 = e~ tuJt <fi(r) , which is a com- 
mon textbook problem on Quantum Mechanics The 
initial scalar profile is then of the form (in dimensionless 
coordinates) 



ipi(0,x < x a ) = -0Q 



0i(O, x > x ) = tpt 



sm(x/a) 
exp 



(11) 



-y/2\U \a 2 -l(x/a) 



The eigenvalues of a are found by imposing the con- 
tinuity of the radial function 4>{r) and its first deriva- 
tive at the discontinuity of the potential, and then are 
given by the solutions to the equation cot(a;o/er) = 
— ^j2\Un\<j 2 — 1, with xq = mRg. Hence, a can be es- 
timated from the argument of the cotangent function as 
(x /(j) ~ (n+l)7r < y/2\Uo\xo, where n = 0, 1, 2, .. labels 
the number of nodes of the initial profile, ipi is evaluated 
from the continuity condition of the radial function at 
x — xq. Finally, the initial mass according to is 



M = 47rm 2 



SVdr^i— (12) 
2 m 
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which becomes an equation that helps to determine V'o for 
fixed n,Ro,M ; a procedure rather convenient since the 
latter are the physical (input) parameters of the system. 

We mention here an additional major simplification 
we can make use of for the SN equations, but not for the 
EKG ones. It is easily seen that Eqs. 0, © an d © are 
invariant under the scaling transformation 

{t, x, U, U 2 , ip} -> [\~ 2 t, A" 1 ^ A 2 f7, A 2 l> 2 , A 2 ^ j (13) 

where A is an arbitrary parameter. A similar scaling 
property was noticed before for the stationary solutions 
of the SN equations, from which one realizes that all sta- 
tionary solutions (for a fixed number of nodes of the field 
are related to each other by a scaling transformation 

By means of (I13|) and an appropriate A, the collapse of 
our fluctuation can be reduced to the study of a con- 
veniently sized configuration concerning /ia£-quantities 
only. Once the hat-configuration has been evolved, we 
apply the inverse transformation to recover the physical 
(no-hat) quantities. 

We focus now on the numerical solution to the SN 
equations. Once the initial profile tpi(0, x) is given, the 
Poisson equation JSJ is integrated with a second order ac- 
curate upwind method inwards under the boundary con- 
dition U(r,Xf) — GM(r,Xf)/xf (being Xf the last point 
of our spherical grid), which is valid at each time slice 
according to the Newton theorems regarding spherical 
objects. The next scalar configuration is determined by 
solving the Schrodinger equation (JJJ using a second order 
finite differences fully implicit Cranck- Nicholson method 
|19| . The procedure is then repeated forward in time. 

To measure the accuracy of our solutions, we use 
the {t, r} component of the Einstein's equations, which 
rephrases the conservation of probability for the Newto- 
nian field ip, 

P := - ^ [rdl(x^) - ^d 2 x {xr)] , (14) 

where (3 is zero for an exact solution. The accuracy 
depends on the time step At and the grid resolution 
Ax, which should be chosen appropriately to assure that 
|A^|/|^| <C 1 in a time step. Thus, we should comply 
with both Ar/(Aa;) 2 < land [At/(Ax) 2 ]\1 + U(Ax) 2 \ < 
1. The former is the condition applied to a free wave 
function, and the latter takes into account the presence 
of a potential in Eq. Q. 

In Fig. [I]we show the runs of an initial 5 (14)-nodc 
scalar fluctuation for a single scalar halo, whose initial 
physical parameters are Mq = 10 11 (1.7 x 10 14 ) Mq and 
Rq = 70 kpc (7.3 Mpc). The corresponding scale param- 
eter is A 2 = 6.38 (0.6) x 10~ 6 , and then a = 14.6 (190). 
The grid spacing is Aa; = 0.25 (4.0) with the bound- 
ary at x f = 1250 (4 x 10 4 ). The time step is Af 
'.') (0.125) x 10~ 2 , and the runs were followed up to a 
physical time of T = 30Gyr. In all cases, \\(3\\ 2 < 10~ 7 . 

We see in Fig. ^that the initial rotation curve profile 
of both evolved configurations is almost flat, but such 
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FIG. 1: Evolution of a 5-node and a 14-node ini- 
tial scalar fluctuation of the form lilt , see the text for 
details. Top: Profiles of the rotation curves v(r) = 
\J GM(r) jr at different times. Bottom: Evolution of the ratio 
Kinetic Energy / 1 Gravitational Potential Energy | (solid line) , 
which shows the virialization of the two systems, and the 
total integrated mass M (dotted line). It is worth to notice 
the so called gravitational cooling after free fall for the first of 
the systems, which is not achieved in the second case within 
the time scale shown, but is expected afterwards. 



flatness is lost at late times during the collapse. The 
smallest configuration, which is supposed to have started 
collapsing earlier within the cosmological picture, rapidly 
virializes and forms a giant Newtonian boson star, though 
not as large as a realistic (dwarf) galaxy. According to 
Eq. i|13|) , the same fate is expected for the largest configu- 
ration, which roughly corresponds to a cluster of galaxies 
fluctuation, but within a time scale much larger than the 
actual age of the universe. 

Some final comments are in turn. The first one is 
about the initial profile JUJ. It is well known that the 
self-gravitating equilibrium configuration of the SN equa- 
tions (the so-called Newtonian boson stars[17| and New- 
tonian oscillatons 0|) are reached by the system after 
the scalar field has had enough time to modify the grav- 
itational potential. We assume that such condition is 
not achieved before turnaround because there are other 
factors involved in determining the scalar profile. There- 
fore, one does not expect the initial profile to be a New- 
tonian oscillaton, but the latter will rather be the final 
configuration of the evolutions shown here. Indeed, Fig.Q] 
confirms the expectations and shows how, in some cases, 
a scalar field fluctuation can virialize and form a New- 
tonian oscillaton. Second, we have seen that a model 
of structure formation within the scalar field hypothesis 
could imply an initial configuration with nodes. The lat- 
ter can be naturally expected since the Compton length 
of the scalar particles is, typically, much smaller than 
the size of the initial configuration, Ac = rn^ 1 <C Rq. 
Then, the wave function can be confined and wrinkled by 
the gravitational well potential (which remains fixed in 
time, see |2jJ, EI) before the non- linear collapse of 
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the fluctuation. 

On the other hand, the scale invariance of the SN equa- 
tions 1|13[) suggests that all galaxies should look similar; 
in particular, some of them could be related just through 
a scaling transformation. This argument can be extended 
to larger and more massive fluctuations, and then it pre- 
dicts that the latter should have similar properties as 
galaxies. 

The SN equations were derived under the assumption 
of spherical symmetry, but we really expect the EKG 
equations to have a tractable non-spherical weak-limit 
too. This would suggest that the full SN equations with- 
out symmetries should be used to explore the realization 
of a realistic non-spherical collapse and to track the col- 
lapse of many scalar fluctuations. Also, the arguments 
about the initial configurations which lead us to can 
be extended to the cases of non-zero angular momentum 
(see for instance pol l22| 1 . 

Another advantage of the SN system is that the scalar 
field galaxy halo can be, in an easy manner, gravitation- 
ally coupled to baryons; the latter would only enter on 
the r.h.s. of the Poisson equation ijSJ 8]. Results on 
this issue and others, like more realistic initial conditions 
and the study of rotation curves profiles, will be reported 
elsewhere. 



Finally, all the results shown here are also valid for 
complex scalar fields, since the dynamics of both type of 
scalar fields is determined by the SN system only, though 
real scalar fields would be preferred from the cosmological 
point of view 0, • 

Summarizing, we have shown how to track the non- 
linear collapse of a matter fluctuation made of real scalar 
fields using the Schrodinger-Newton system, and the kind 
of properties it should have before and after separating 
from the cosmological expansion. 
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